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Abstract 

We consider a model of nonlinear superfluid flow past a periodic array of 
point-like scatterers in one dimension. An application of this model is the 
determination of the critical current of a Josephson array in a regime appro- 
priate to a Ginzburg-Landau formulation. Here, the array consists of short 
normal-metal regions, in the presence of a Hartree electron-electron interac- 
tion, and embedded within a one-dimensional superconducting wire near its 
critical temperature, T c . We predict the critical current to depend linearly 
as A(T C — T), while the coefficient A depends sensitively on the sizes of the 
superconducting and normal-metal regions and the strength and sign of the 
Hartree interaction. In the case of an attractive interaction, we find a fur- 
ther feature: the critical current vanishes linearly at some temperature T* 
less than T c , as well as at T c itself. We rule out a simple explanation for the 
zero value of the critical current, at this temperature T*, in terms of order 
parameter fluctuations at low frequencies. 



1 



Typeset using REVTeX 



I. INTRODUCTION 



The problem of nonlinear flow, as modeled by a nonlinear Schrodinger equation in the 
presence of a scattering potential, relates to a variety of physical situations. One example is 
a weakly interacting Bose gas with impurities, while another is a Josephson array in a regime 
appropriate to a Ginzburg-Landau formulation. Here, the array consists of short, normal- 
metal (N) regions in the presence of a Hartree electron-electron interaction, and embedded 
within a superconducting wire near its critical temperature, T c . The latter example falls 
within the rapidly developing field of study of the interplay of the proximity effect and 
charging effects due to interactions in disordered, inhomogeneous superconducting systems 
|jl|-0|. Further examples arise within nonlinear optics and in the study of gravity waves 
on deep water [|7J. 

We will concern ourselves with the existence of time-independent solutions of the flow 
equations in the presence of a supercurrent. In general the value of the supercurrent is 
limited by a maximum value, the critical current, above which the flow becomes unsteady 
and a time-dependent solution must be sought. Recently Hakim || considered the superfluid 
flow past a single, repulsive scatterer in one dimension and deduced the dependence of the 
critical current on the scattering strength. We extend this model to the geometry of a 
regular array of scatterers. We find a rich dependence of the critical current on the scatterer 
separation and strength, with markedly different behaviour in the cases of repulsive and 
attractive scatterers. 

Although we will focus here on the steady-state flow, the physics of unsteady flows in 
related problems has also received attention. For instance, unsteady flow past a single 
scatterer in one dimension || represents a transition state towards the emission of solitons 
and hence phase-slip nucleation. It has also been pointed out |§ that steady-state flow 
may not be stable to mechanisms of quantum tunnelling, leading to phase-slip nucleation 
and non-zero dissipation even below the (quasiclassical) critical current. Higher dimensional 
analogues of the problem are also of theoretical and experimental importance, although less 
tractable analytically: recent work has simulated numerically the flow of a superfluid past an 
obstacle |TU[ or through a constriction |TT|] in two dimensions, for which vortices are nucleated 
above the critical supercurrent. These investigations relate directly to experimental work 
[ T2| on superfluid 4 He. 

As mentioned above, an important application of this model is in the description of su- 
percurrent flow in a mesoscopic SNS device near the critical temperature of the S region. The 
Ginzburg-Landau equations for the order parameter near the critical temperature of such a 
device are equivalent to the flow equations for a superfluid in the Hartree approximation. 
Furthermore, a short N region may be modelled as a point-like scatterer whose strength 
diverges as T c is approached and the Ginzburg-Landau correlation length diverges. We will 
exploit this equivalence to translate our results directly into experimental predictions for 
the temperature dependence of the critical current of a Josephson array, consisting of short 
N regions embedded in a superconducting wire near T c . 

The dependence of the critical current on the temperature is a property that has already 
been examined in a number of related systems. For example, it is an established result 
|rf3| that a superconducting wire (with no N regions) has a critical current that behaves as 
(T c — T) 3 / 2 as T c is approached, while the introduction of a single, short N region leads to a 
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critical current that depends quadratically as (T c — T) 2 [p]]. In addition it has been possible 
to fabricate and measure experimentally a superconducting 'microladder' [14|, consisting of 
a pair of S wires with connecting side-branches: in this geometry, a well-defined correction 
to the (T c — T) 3 / 2 dependance of the critical current of the superconducting wire has been 
observed and explained theoretically. Courtois et al. [15j have also measured the critical 
current for another geometry (somewhat complementary to that considered by us) of an 
array of short S regions near T c within a N metal, again with a reasonable fit to the theory 

The case of two short N regions, comprising a double-barrier structure, has also been 
examined theoretically, by Zapata and Sols . They interpret such a system as representing 
a nonlinear analogue of resonant tunnelling: such a scattering structure is well known in a 
linear system to lead to phenonema such as the sensitivity of the transmission coefficient to 
the scatterer separation, peaking near well-defined resonances due to multiple inner reflec- 
tion. The introduction of nonlinearity leads immediately to markedly different behaviour, 
the added complexity precluding even any kind of crossover regime in terms of a nonlinearity 
parameter. The critical current of the double-barrier structure depends as (T' c — T) 1 / 2 , where 
T' c is some temperature below T c . 

For the geometry of an array of short N regions embedded within a superconducting 
wire, we find a linear dependence of the critical current as A(T C — T), in contrast to the 
cases described above. Furthermore, the coefficient A depends sensitively on the sizes of the 
S and N regions and the strength and sign of the electron interaction. The calculation of the 
chemical potential, while trivial in the geometry of a finite number of scatterers, becomes 
more complex in the array geometry, and follows by the incorporation of a normalization 
condition to fix the total boson number. This procedure leads to the increased complexity 
of the flow solutions. 

We find a further feature in the temperature dependence of the critical current in the 
case of an attractive interaction: the critical current vanishes linearly at some temperature 
T* less than T c , as well as at T c itself. We will rule out a simple explanation for the zero 
value of the critical current at this temperature, T*, in terms of order-parameter fluctuations 
at vanishingly low frequencies. 

The plan of this paper is as follows. In Section [□], we describe the model and flow 
equations for a superfluid in the presence of scatterers. In Section [I il| we review briefly 
the solution of these equations for a single scatterer in the time-independent regime. In 
Section |TV] we address the geometry of an array of scatterers and present the results for 
the critical current as a function of scatterer strength and separation. In Section [V] we 
derive analytically the form of the critical current in the various limiting cases of scatterer 
strength and separation, while in Section [V| we translate these results into predictions for 
the temperature dependence of a Josephson array near T c . Section |VII| includes a discussion 
on order-parameter fluctuations at the special separation value of L = 2\g\, at which the 
critical current vanishes, and Section |VIII| concludes with a summary. 
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II. MODEL 



We derive the relevant equations first for the problem of a one-dimensional superfluid 
in the presence of delta-function scatterers. Taking the scatterers to have strength g a and 
positions r a , the Hamiltonian may be written 

H = E (-^ d l+J29J(x i -r a ) \ + -x k ), (1) 

i \ ot I i^k 

where i labels the N bosons, which interact via a short-ranged, pairwise potential of strength 
A. We employ the Hartree approximation W?\ to write the ground state wavefunction in the 
symmetrised form, 

1 N 

X 2 , ...,X N ,t) = — ^ Yi ^fai, *)> 
L,sys i=l 

where C sys is the system size. This leads to the following nonlinear Schrodinger equation 

for i(>(x, t): 

~ Wjl> + M V + E 9J(x - r a )^ = idt^ (2) 

1 a 

where we have rescaled length according to units of the healing length, 1^ = h/(Xnm) 1 ^ 2 , 
and energy (and the g a ) according to units of the Hartree energy, nX, where n is the average 
particle density. The speed of sound in the condensate in the absence of the scatterers, ie. 
(nA/m) 1//2 , is equal to unity with this choice of units. 

In addition, we must enforce the condition of fixed total boson number, which is achieved 
by normalizing the wavefunction according to 



{x)d,X — C-sys- (3) 

We now introduce the number-phase representation (the Madelung transformation [|T8|j): 

i/>(x,t) = y/p(x,t)exp(iS(x,t)). (4) 

Inserting eqn. (H) into the Schrodinger equation, eqn. (j^), gives the continuity equation and 
the Bernouilli equation respectively, 

d t p + d x (pd x S)=0, (5) 
+ dtS + +Y,9aS(x -r a )+p = 0, (6) 

while the normalization condition, (§), becomes 

p(x)dx = £ sys . (7) 
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For time-independent flow, we may set dtp to zero, while the phase S(t) advances uniformly 
in time according to the Josephson relation, d t S = — p, where p is the chemical potential. 
The continuity equation, eqn. (Q), then integrates to 

pd x S = j, (8) 

where j, a constant, is the supercurrent density Writing <p(x) = p(x) 1 / 2 , the flow equation, 
eqn. (BI), becomes 



\& x <i> + 3 - \4 + 9aS(x - r a )<P + il = 0. (9) 
2 ^ 20 d 



We will discuss the solution of this equation at length in the following sections. We see that 
it is of a general Ginzburg-Landau form, as we clarify further in section [V]]. We consider 
first the case of a single scatterer before turning to an array of scatterers. 



III. SINGLE SCATTERER 

In this section we review the results for the critical current of a single scatterer |§. We 
take the scatterer to be of strength g and placed at a position x = 0. The flow eqn. (||) gives 
the jump condition 

\[d^\t=gm- (io) 

In addition, from the normalization condition, eqn. (0), we enforce that 

4>(x) — > 1 as x — > oo. (11) 

From eqn. (|TT|) and the flow equation, eqn. (H), we find the chemical potential immediately 

as 

With the above relations, the flow equation may be integrated as follows: 

p(x) = f + (i - j 2 ) { ^ j [vT=rW + «1 , { ! ^ °' (12) 



<?<0. 



The integration constant, a, still needs to be determined from the jump condition, (10). 
For an attractive impurity, g < 0, the jump condition may be fulfilled for all values of the 
supercurrent up to the speed of sound, j < 1. This is connected to the fact that the value of 
p at the impurity (p(0)) may become arbitrarily large. Thus the critical current for a single 
impurity, j C) o, is in this case 

jc,o(9 < 0) = 1. 
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FIG. 1. Critical current, j c , for a single, repulsive impurity of strength g > 0. 



In contrast, for a repulsive impurity, g > 0, the jump condition can no longer be satisfied 
for all j up to 1, due to the restriction that < p(0) < 1. Instead we have that j c ,o{g > 
0) < 1. Although Hakim || obtained an implicit formula for the critical current, we obtain 
in Appendix [A] a relatively simple explicit formula. The resulting dependence of j c ,o{g) on 
g is illustrated in figure |1[ In the limits of large and small scattering strength, g, we have 

Jcfliff > 0) - j ! _ |( 2 ^)i/3 ) g i (13) 

Note that the result for large g will find relevance in section [Vl| when we examine the 
equivalent problem of the critical current in an SNS junction (with a single normal-metal 
region) close to the critical temperature. Before elaborating on this interpretation, we keep 
to the example of a superfluid, and turn to the case of a periodic array of scatterers. 



IV. ARRAY OF SCATTERERS 

We now examine the case of a periodic array of scatterers, of equal strength g and 
separation L. We place them at positions r a = {2a + l)L/2, so that <fi(x) is symmetric 
about x = and we may restrict attention to \x\ < L/2. We first show how the flow 
equations may be integrated before discussing the results of their full solution. 

A. Integration of the Flow Equations 

The first integral of the flow equation, (||), may be written as follows: 

{d^f = 4 - 2^ 2 0(O) 4 + 2/i0(O) 2 + r 



m 2 

We may now factorize the right-hand side of this expression: this operation greatly simplifies 
the following analysis. We find 
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(d x <j>) 2 = -W 2 - 0(O) 2 )(0 2 - «)(0 2 - 0), (14) 



where 



1/2 

p(o) , ^ - ! 

- /i — ± 



p(0) 



(15) 



This leads to the following solution for p(x): 



p(x)=P+(p(0)-P)\ S f )(V^Px + ReK(k),k), 



fc 2 sn 2 



where the symmetry of p(x) about the origin has been automatically incorporated, sn and 
K are the elliptic integral and complete elliptic functions respectively |19[], and 



It still remains to determine the two integration constants, p and p(0). These are specified 
by the normalization condition corresponding to eqn. (|7|), that is, 

2 r L / 2 

p(x)dx = 1, (17) 



LJo 

and the jump condition, 



The latter condition may be rewritten, using the first integral in eqn. ([14]) , as follows: 

?TO^ 2 ?(l7ra + 2 ' 1 -' )(i/2) -' )(0) ' (19) 

This represents an implicit equation for j, since the right-hand side is dependent on j 
through p and p(0). The problem of calculating the critical current of the array has now 
been reduced to finding the maximum value of j for which the two conditions fllTD and ( |j"U| ) 
may be satisfied simultaneously. 



B. Results 

Having described the integration of the flow equations, we present here the results for 
the critical current of the array. Further details of the working towards these results will be 
described in the following section. 

It is clear that in the limit of L — ■> oo, j c must approach its single-impurity value, j Cj o: 

jc^jcM, L -»• oo. (20) 
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FIG. 2. Critical current, j c 
and (b) an attractive array (g = 
dotted line) for large L. 



(solid line), as a function of L for (a) a repulsive array (g = 5) 
- — 1.0). In each case, j c approaches the single- impurity value (the 



For general values of the parameters g and L, the determination of j c must be performed 
numerically. Figure |2| shows j c as a function of L for a typical repulsive and attractive case. 
We see that the correct behaviour is reproduced in the limit of L — > oo. We also see that 
in the attractive case, the critical current vanishes altogether at one special value of L. We 
will show below that this value is given by L = 2\g\. We also see that j c diverges at small 
L in both the repulsive and attractive case. 

In the limit of small scatterer separation, L — > 0, we find further that j c satisfies the 
following scaling form: 

j c ^Lf(gL), L^0. (21) 

The scaling function / possesses a surprisingly nontrivial structure which we have obtained 
numerically and show in figure |^. Notice the marked difference in these forms for the 
repulsive and attractive cases. We see also that, with interpretation of this system as a 
superfluid as described in section |T|, the scaling variable \g\L is unaffected by the energy 
and length rescaling (in units of the Hartree energy and healing length respectively) and 
hence of the boson interaction A. 

The limits of large and small separation needs a more precise definition, which in fact 
differs in the repulsive and attractive cases: the limits correspond to L ^> L and L ^ L 
respectively, where L = 1 for the repulsive case and Lq = \g\ for the attractive case. 

The scaling function / will find direct relevance in the prediction for the temperature- 
dependence of the critical current of a Josephson array in section |VJ In the following section, 
we will demonstrate analytically the following limiting behaviour of this function: 
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The critical current j c for small L then follows by s = \g\L and the scaling relation (^T|). We 
will also show analytically that the critical current vanishes as the special point L = 2\g\ is 
approached, in a linear fashion as j c oc |L/2|p| — 1|. 

V. LIMITING FORMS OF THE SOLUTIONS 

In this section we derive analytically the limiting forms of the scaling function given 
by eqns. (|22| ) and (0), and hence of the critical current at small separations. We take 
the repulsive case first, which requires the limit L < 1, and then the attractive case, 
which requires L <C \g\. Note that the condition itself for small separation is different in 
the repulsive and attractive cases. While staying within this condition, we will examine 
separately the subcases of \g\L <C 1 and \g\L 3> 1. 

In general, we find at such small separations the chemical potential approaches a large 
positive (negative) value in the repulsive (attractive) case, reflecting the large 'potential 
energy' of the scatterers. In addition, for small \g\L, the density p(x) remains close to 1 for 
all values of x. In contrast, for large \g\L, the value of p(x) approaches 2 at the origin and 
at the scatterers for the repulsive case, and at the origin and \g\L ^> 1 at the scatterers 
in the attractive case. 

To derive the critical current, it is necessary to identify the somewhat subtle interplay 
of the various parameters of the problem, which requires a separate and quite different 
discussion for each of the four cases. We also include a discussion of the special points 
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L = 2\g\, at which the critical current vanishes entirely. Our approach in each case will 
be to indicate the existence of a local (rather than global) maximum in the supercurrent. 
Strictly speaking, we still need to justify these choices of a local maxima as the relevant 
values for the critical current. To do so, we appeal to the numerical results of section 
IV B| which establish the smooth interpolation from the limit of small separation to large 



separation and hence to the single scatterer result for the critical current, which certainly 
does represent a global maximum. 



A. Repulsive case at small separation (L <C 1) 

In the repulsive case, the limit of small separation requires that L < 1. We examine 
separately the subcases of gL <C 1 and gL ^> 1, while staying within the limit of L <C 1. 

In both subcases, we have that /i> 1, and hence k <C 1 by eqn. (|16|). The sn function 
is then well approximated by a standard cosine function: 



p(x) =(3 + (p(0) - (3) cos 2 (J2px) 



(24) 



where (3 = j 2 /(2pp(0)). We also have that p(x) must contain no more than one half- 
oscillation between impurities, so that ^/2p,L/2 e (0,7r/2). The normalization condition, 
eqn. (|17]) consequently simplifies to the form, 



(3) 



1 + 



2pL 



■ sin 



2pL) 



(25) 



As shown in Appendix [B], this limiting form leads to a critical current which behaves as 
j c ~ tt/L for gL <C 1 and j c ~ n 2 / (gL 2 ) for gL 3> 1 (while 



B. Attractive case at small separation (L <C \g\) 

In the attractive case, the limit of small separation requires that L <C \g\. Note that this 
condition itself is separate from that in the repulsive case. Again, we will concern ourselves 
with the two subcases \g\L <C 1 and \g\L ^> 1, while staying within L <C \g\. In both 
subcases, the chemical potential is large and negative: \p\ 3> 1, p < 0. The behaviour of k 
however is different according to the limit of \g\L: we have k ^> 1 for \g\L 1 but k ~ 1 
for \g\L ^> 1. The sn function then reduces to either a cosine function or a tanh function: 



p(x) = S /3 + (p(0)-/5)sec 2 ( v /p(0)-/5x), |(?|L«1, m 
\f3+(p(Q)-f3)coth 2 (p-TS^px), |^|L>1, 

P=hog(-A-r)- (27) 



2 °\ife-l 



Further details of the derivation of the critical current are contained in Appendix 0: we find 
j c ~ (2|(?|/L) 1 / 2 for \g\L 1 and j c ~ 2g 2 L exp(— \g\L) for \g\L ^> 1. This concludes the 
derivation of the behaviour of the critical current and hence the scaling function contained 
in eqns. (|22|) and (JEJ. 
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C. Attractive case in the limit L — > 2|<?| 



In this section, we consider the attractive case at separations near the special value, 
L = 2\g\, at which the critical current vanishes entirely. This limit is characterized by a 
divergence in one of the parameters, namely k — > oo. Consequently, p — > p(0)/2, and hence 



p(x) - p(0)sec 2 (^(0)x) 



L — > 2|g|. 



(28) 



The normalization condition, eqn. (0), and the jump condition, eqn. flT9]) , become equivalent 
in this limit, and simplify to: 



L 

2 



p(0) tan 



P(0) 



L 



(29) 



Eqn. (|29| ) for p(0) is soluble for all values of L. 

At this point it is still not obvious that the special point L = 2\g\ does not admit a 
steady solution at any non-zero supercurrent. To demonstrate this fact, we consider the 
point L = 2\g\(l + e) for some small e (positive or negative), and show that j c vanishes as 
|e| -> 0. 

In this limit, we identify the small parameter 7 = p — p(0)/2, I7I <C 1. We will find 7 to 
be proportional to e at the critical current. Then 



a 



7± 7 2 " ' 



1/2 



P(0), 



We find that k 2 = p(0)/(2( / y — f3)), and hence the limiting form for the density is 



p(x) =(3 + (p(0) - (3)sec 2 



p(0)-/3 1 



2p(0) 



x 



(30) 



As we show in Appendix 0, this limiting form leads to a critical current j c that vanishes 
linearly in |e| as claimed. 

In section |VII| we will examine order parameter fluctuations, propagating as sound waves, 
at this separation L = 2\g\. We will rule out the existence of such fluctuations at arbitrarily 
low frequencies as a simple explanation for the zero value of the critical current at this value 
of the separation. Instead we may view this effect as a consequence of the nonlinearity of 
the flow equations. 



VI. JOSEPHSON ARRAY NEAR T c 

Having the determined the behaviour of the critical current for the periodic array of 
scatterers, we now show how these results may be translated directly in experimental pre- 
dictions for the temperature dependence of the critical current of a Josephson array in a 
regime where a Ginzburg-Landau formulation is appropriate. Here, the array is near the 
critical temperature, T c , of its S regions, and has the quasi-lD geometry shown in figure £|: 
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N 

FIG. 4. The geometry of the Josephson array. 



the S regions are of length L$ while the N regions are of length <C L$ and subject to a 
Hartree potential, V. 

We first show that the Ginzburg-Landau equations for the array are of the same form 
as the flow equation, eqn. (||), for a superfluid. Furthermore, as shown by Zapata et al. @, 
the short normal-metal regions, in the presence of the Hartee potential, may be mimicked 
by a series of repulsive or attractive delta-function scatterers, depending on the sign of the 
interaction V. 

The Ginzburg-Landau equations for the order parameter \fr(x) of the S region may be 
written as Iffi! 

h 2 

-—dim + a(T)m + 6|^| 2 ^ = 0, 
4m 

where a(T) is negative and proportional to (T c — T), b = 2TCK 2 e 2 h 2 /m 2 c 2 and k is the 
Ginzburg-Landau parameter. Writing \1> = |^| exp(ziS') and defining the supercurrent (per 
unit of cross-sectional area) as j = (eh/m)\ip\ 2 d x S, we have 

-^ + «W + ^ 3 + £^ = 0. (31) 

In addition to the above differential equation for ip, we need to specify appropriate boundary 
conditions. In the absence of any N regions, the order parameter of the homogeneous 
superconductor, and hence condensate density, is fixed at the value ipsiT). Furthermore, 
for a geometry containing only a single N region, the appropriate boundary condition is to 
set that ip{x) approaches ips f° r limiting values of x far from the N region, by analogy with 
the boundary condition used in section fTD] in considering a single scatterer. For an array, 
however, no such condition at infinity can apply; instead, we fix the total boson number at a 
certain temperature to equal that in the absence of the N regions. This gives the condition 

1 /-is/2 

— / ip(x)dx = ip s , 

L s J-Ls/2 

in similarity to the normalization condition, eqn. (|TTD, of section [IV]. Here the integration 
extends over a single S region. To proceed, we rescale the position variable and current as 
follows: 
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167TK 2 e£(T) 
where the coherence length is defined as 

n 



^ h ° 2 T 



£(T) 



2(bmy/ 2 ip s 
/ T\~ 1/2 

a I 1 " ^ 



and = 0.74^o for a clean superconductor and £q = 0.85v^£o f° r a dirty superconductor, 
where £o and £ are the zero-temperature coherence length and mean free path respectively. 
In these rescaled variables, the Ginzburg-Landau equation, eqn. (|31~1) , becomes 



1 T 2 
2 



-^0-^ + 3 + — = 0, (34) 



where /x = |a|/(6 , ?/ ; ^o)- We see that for the geometry of a single N junction, the chemical 
potential becomes \i = 1 + J 2 /2. We also see that eqn. ([34]) reproduces the flow equation 
(|) as expected, together with the normalization condition (|T7j). 

Notice that the rescaling leads immediately, by eqn. (^), to a (T c — T) 3 / 2 dependence 
for the critical current of a S wire with no N regions, in accordance with established theory 



see e.g. [0]). If we now include one or more short N regions, the flow equation in an N 



region in the rescaled variables, corresponding to eqn. fl3~4j), reads as 

1 „, . „ . J 2 



where 

e = v/(hp 2 s ) 



2 



= sgn(^)^, 

and £jy = U 2 / (Am\V\) is the Ginzburg-Landau correlation length in the normal-metal. Fol- 
lowing Ref . , we may now integrate trivially the flow equation fl3"o] ) over the whole N region 
(located in y G (y~,y + ), say) to give 

-\^-=^f^ (36) 



as long as 
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Notice that we have dropped the current- dependent term on the right-hand side of eqn. (p5|) 
as is consistent for a sufficiently short N region near T c (see 0). By identification of eqn. (|36D 
with eqn. (|i~8|), we see that we may model the short N regions by delta-function scatterers 
of strength 

g = sga{V) . (37) 

At the same time, following the rescaling of the size of the S regions, Ls, under the trans- 
formation fl32|) , we may also identify the scatterer separation as 

L = ^ . (38) 

By eqns. ( |3"T| ) and Q3HD, we see that the limit of T — > T c , and hence £(T) — > oo, corresponds 
to taking g — > oo and L — > simultaneously, such that the product pL(= L N L S / {2^)) is 
fixed. 

Given the above information we are now in a position to predict the temperature depen- 
dence of the critical current of the array. For orientation however we start with the simpler 
geometry of a single S-N-S junction with a repulsive interaction. As the temperature ap- 
proaches T c , g diverges by eqn. ([37]) and we insert the large g limit of the critical current 
for a single impurity, eqn. flT3]): J c = \/{2g). Relation (|33|) then gives the critical current of 
the SNS junction as 

. me N f t 

per unit area, ie. it varies quadratically as (T c — T) 2 . 

We turn now to the array. As the temperature approaches T c , we have L — > and we 
become able to apply the scaling form (|21~D for the critical current per unit area: 



3c(T) = h ° f{s) (l-^r 
LnLs 



2a 



Recall that f(s) is the scaling function discussed in section [TV], illustrated in figure |3| for the 
repulsive and attractive cases and obeying the limiting forms contained in eqns. ( |22| ) and 
(p3|). We see that we have a linear dependence of the critical current as j c = A(T C — T). In 
addition, the associated coefficient, A, of this linear dependence depends on the parameters 
L N , L s and V through the product L N L S V: 

he 2 / 2mL N L s V 

Notice that both large and small values of this product may be probed while staying within 
the requirement that L N <^ £ N . 
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In moving away from T c for the attractive (V < 0) case, an addition feature arises in 
accordance with the vanishing of the critical current at L = 2\g\: we have that the critical 
current vanishes at the temperature T* < T c , as well as at T c itself. Here, T* is determined 



Moreover, since j c vanishes linearly in |L/2|gr| — 1|, we have that the critical current vanishes 
linearly in \T — T*\ as T* is approached. 

This concludes our discussion of the temperature dependence of the critical current of 
the Josephson array near T c (and T*). We now focus on the case of T — > T* to examine 
whether order parameter fluctuations exist at arbitrarily low frequencies in this limit. 



Having identified the special value (L = 2\g\) of the (attractive) scatterer separation at 
which the critical current takes a zero value, we here examine whether a simple explanation 
for this effect exists in terms of order parameter fluctuations at arbitrarily low frequencies. 

We use a standard procedure to describe such fluctuations in a condensate of non-uniform 
density (see, for example, Giorgini et al. f21fl ): we perturb the wavefunction ip(x,t) (see 
Section [Tl]) by a small contribution that is oscillatory in time: 



The flow equation @ may then be linearised in these small oscillations to give the following 
coupled flow equations for u(x) and v(x): 



by 




VII. ORDER PARAMETER FLUCTUATIONS AT L = 2\G\ 



ij,(x, t) = e-^ <f){x) + u{x)e~ iwt + v*{x)e 



Cu(x) + p(x)v(x) = uu(x), 
p(x)u(x) + Cv(x) = —uv(x), 



(39) 
(40) 



where 



£ = -d 2 x /2 -p + 2p(x) + y £g a 6{x - x a ). 



Q 



These equations are supplemented with the normalization condition, 




(41) 



while trivial integration of the flow equations (39) and flUf) over each scatterer gives the 



jump conditions, 




(42) 
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Our aim is to determine whether fluctuations in the wavefunction, propagating as sound 
waves, exist at arbitrarily low frequencies at the separation L = 2\g\. To this end, we search 
for a consistent solution for u(x) and v(x) in the limit of uj — > 0. 

For preparation we review the solutions in the homogeneous case, ie. in the absence of 
any scatterers. In this case, the perturbations are plane-waves, u(x) = exp(—ikx)u/ C sys and 
v(x) = exp(— ikx)v / £ sys , where £ sys is the system size as before. The flow equations (|39| 
and pop, together with the normalization condition, fl¥ip, are now easily solved to give 



u 



C + uj 

2uj '' 
C — UJ 

2uj 



where C = k 2 /2 + p. The dispersion relation reads 

1 



(43) 
(44) 



UJ 



-k 2 (k 2 + 2p). 



We now generalise to the case of the non-uniform density at the separation, L = 2\g\. 
In this case, we need to incorporate the jump condition given by fl4"2| ) at each scatterer. 
Given the periodicity of the arrangement of the scatterers, we will search for purely periodic 
solutions for u and v. 




x + L) 




(45) 



so that we may restrict attention to only a single region \x\ < L/2 with one jump condition. 
Eqn. fl4"5|) may be seen as a Bloch theorem for u and v, at Bloch momenta which are precisely 
multiples of n/L: odd and even multiples give rise to solutions in u and v which are odd 
and even in x, respectively. 

Bearing in mind the form of the solutions (f43D and (44) for the uniform case, a consistent 
series expansion for u(x) and v (x) in the limit of uj — ► becomes: 



u\x) 



U {X) 



;i+0(uO), 



The normalization condition, (|4l|) , to zeroth order then gives Mq = Vq\ by analogy with 
the homogeneous case, we take u = v . The zeroth-order differential equation for u (x), 
corresponding to eqn. (|39D , becomes 



U[X) 



for x E (—L/2, L/2). We substitute the form ( p8|) for the density p(x): we find 



6 sec 2 ?/) u(y) = 0, 



(46) 
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where y = J p(0)x. The two independent solutions of this differential equation are as follows: 



sm y 
cos 2 y ' 



n 2 (!j) = /•' ( 1,2; ^;cos 2 y) cos 3 ?/, 



where F is Gauss's hypergeometric function |19fl . By the Bloch theorem, eqn. (|45j) , we will 
take u{y) to be either purely odd or purely even: u{y) = Au\{y) or u{y) = -Bw 2 (y). We now 
check whether either solution is compatible with the jump condition. We see that this check 
is independent of the coefficients A and B. Taking the odd solution first, Ui(y), the jump 
condition ([42]) becomes 

L 
2 



^Jp(G) [cot(y ) +2tan(i/ )] 



where yo = Jp(0)L/2. This simplifies by use of eqn. fl29|) to 



= p(0) + 



which clearly cannot be satisfied for any L, as the right-hand side always exceeds zero. This 
solution must therefore be discarded. We are left with the even solution, U2{y), for which 
the jump condition fl42|) becomes 



5cos 2 (y ) 
sia(y )u(y ) 



cot(y ) + 2tan(?/o 



which follows from standard properties of the hypergeometric function [19]. Use of eqn. (29) 
simplifies this condition to the form 



f(i,2;^cos 2 (i/o)) =5. 



(47) 



However it is easy to check that < F(l, 2; 7/2; z) < 5 for < z < 1, with F(l, 2; 7/2; 0) = 
and F(l,2;7/2;1) = 5: clearly the condition (|47| ) cannot be satisfied either for any L > 0. 

We conclude that a periodic solution for u(x) does not exist in the limit of u> — ► and 
hence such order parameter fluctuations do not exist at arbitrarily low frequencies. Instead 
it seems we must appeal to the nonlinearity of the flow equations as an underlying cause for 
the surprising effect of a zero critical current at L = 2\g\. 



VIII. SUMMARY 

In this paper we have examined the superfluid flow past an array of point-like scatterers 
in one dimension. We have determined the critical current of the flow, above which the 
flow becomes unsteady. While the result for a single scatterer is recovered in the limit of 
large scatterer separation, we find a scaling form for the critical current in the opposite 
limit of small scatterer separation. The scaling function takes a particular form, separate 
in the repulsive and attractive cases, that we have obtained numerically, as well as derived 
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analytically in the various limiting cases of scatterer strength and separation. We also find 
the additional feature in the attractive case that the critical current vanishes altogether at 
one special value of the scatterer separation (L = 2\g\). 

While these results are applicable to a variety of physical situations, an important ap- 
plication is in the prediction of the temperature dependence of a Josepshon array, in the 
presence of a Hartee potential and near T c . In contrast to dependencies already derived and 
observed experimentally in other geometries, we find for the array a linear dependence of 
the critical current as A(T C — T). The coefficient A depends sensitively on the size of the 
normal regions (L N ) and of the superconducting regions (L s ) and the Hartree interaction, 
V, through the product L S L N V. In addition, for the attractive case (V < 0), the critical 
current is suppressed to zero as T — * T*, as well as at T c itself, where T* is some tempera- 
ture less than T c . We have ruled out a simple explanation for this suppression of the critical 
current in terms of order parameter fluctuations at low frequencies, instead appealing to the 
nonlinearity of the flow equations as an underlying cause. 
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APPENDIX A: SINGLE IMPURITY: CRITICAL CURRENT 



To find the critical current for a single impurity, we replace the integration constant 
a in eqn. ( |L2"D in favour of the constant p(0), where tanh 2 a = (p(0) — j 2 )/(l — j 2 ), and 



reformulate the jump condition ([II]) as 

^ 2 p(0) 2 = (p(0)-j 2 )(p(0) 



1; 



where we have used the first integral of the flow equation, eqn. 
j Cj o, is determined by the condition dj c /dp(0) = 0, or 

2/p(0) = (p(0)-l)(3p(0) 



1 



2jen) 



Eliminating g from eqns. 



and (I 



we find 



(Al) 

The critical current, 
(A2) 



P(0) = ~ [-1 + (1 + 8j 2 ) 



(A3) 
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while eliminating j in a similar manner gives 

(p(0)-l) 3 + 2<? 2 p(0) = 0. (A4) 
Eqns. ( |A3| ) and ( |A4j ) give the following solution for j c : 



, 2 2# 2 3# / i2 1 \ 2 /# 2 1 N 

Jo,o(g >o) 2 = i- -f + 4 -- + -+/ - + - 



R 3 



3 V 3 RJ * V 9 ,R 2 / 

2^ 



-1 + ^1 + 8^/27 



This formula reduces to the forms shown in the main text in the limits of large and small 
scattering strength, g. 

APPENDIX B: REPULSIVE IMPURITIES AT SMALL SEPARATION 

In this section we obtain the critical current for repulsive impurities at small separation, 
Lcl. We treat the subcases of gL ^> 1 and gL separately. 

1. The subcase gL^> 1 

In this limit, the have that the density p(x) approaches 2 at the origin and zero at the 
impurities. Using eqn. (|24"D, we see that the latter limit means that y/2jlL/2 ~ n/2. Writing 



%4 = !-7, (Bl) 



where 7 < 1, we have from eqn. (£4|) that 



p(L/2) ~ ^ + 2 7 2 . (B2) 



2tt 

In addition, the normalization condition, eqn. (fffi) , leads to p(0) = 2(1 — 27/71*). Inserting 
eqn. ( |B2"D into the jump condition, eqn. ((T^), we find 

j 2 ^-/p(L/2) 2 + 4pp(L/2) 



2tt 2 " J L 2 

or 

'2 



^(^ + 2 7 ) 2 +j 2 + 
2^7 



7T 



^L 2 



1-1- 



7T 



21 1/2 



Maximisation of this expression with respect to 7 is now trivial and gives 

7T 2 

3c = -£2, ^ > 1, 

as required. At this critical current, we have 7 = tt/ (2gL) <C 1 and p(L/2) = n 2 / (gL) 2 <C 1. 
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2. The subcase gL <C 1 



We now take the opposite limit of gL <C 1. In this case, we have that the density p(x) 
is close to 1 for all values of x. We set u = p(0) -1 < 1, and define 7 as before, in 
eqn. (|BT|). As may be verified at the end of the calculation, we have that «,7 < 1 (for j 
near its critical value). We first use the normalization condition, eqn. (|25|), and the jump 
condition, eqn. flT§|), to determine a relation between u and 7. This will allow us to obtain 
an expression for j in terms of only 7, which may then be maximised simply. 

The normalization condition, eqn. fl25|), now reads 



1=P + 



1 



2 7 



7T — 27 



leading to 



u + ... 



Since /3 = j 2 /(2pp(0)), this gives 

j 2 L 2 



7T Z 



1_ 12 

7T 



« Z + . . . 



(B3) 



(B4) 



In eqn. ( |B~3| ) and ( |B4| ) we have kept to first order in 7 and second order in u, which may be 
checked to be an appropriate level of accuracy at the end of the calculation. The limiting 
form for the density, eqn. (23) now gives 

p(L/2) =(3+ (p(0) - (3h 2 
+ 2«7 2 . 



Inserting this into the jump condition, 



J 



2pp(0) 
0), we find 

f-+j 2 + 4/iM7 2 . 



Notice the cancellation of the j 2 terms, which allows us to find 



1 2tt 



(B5) 



Inserting eqn. ( P5|) into eqn. (|B4|), we find 



fL 2 



1 - 



4 7 (gL)' 



+ 



H'- 7T 47T 2 7 2 

This expression may now be maximised simply with respect to 7, with the result 

2/3- 



3c 



71 
I 



2 U 2 / 



(B6) 



as required. At the critical current, we have that 

1 



7 



27TV3 



(gLf 3 « 1; 



u 



IT 



2/3 



(gL) 1 " « 1. 
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APPENDIX C: ATTRACTIVE IMPURITIES AT SMALL SEPARATION 



In this section we obtain the critical current for attractive impurities at small separation, 
L -C \g\- We again treat the subcases of \g\L 1 and \g\L <C 1 separately. 



1. The subcase \g\L ^> 1 



In the limit of \g\L 3> 1, we have that p(x) is close to zero at the origin and p(L/2) ^> 1 
at the impurities. As may be confirmed at the end of the calculation, we have \/3\ ^> \a\ and 
hence k ~ 1, or more precisely, 



where a = — j 2 /(2pp(0)). The jump condition (|TP|) gives 

0~<7 2 -2|p|, 

and hence p(x) = g 2 csch 2 (p — The normalization condition, eqn. (|T7l), then leads to 

l = lMexp(-2p+| S |L), 



and hence p(L/2) = \g\L ^> 1. The expression for k, eqn. (|CT|), together with the definition 
of p, eqn. fl2"T|), together give 

f = g 2 p(0)(4\g\LeM-\g\L) - p(0)), 

which is again easily maximised with respect to the parameter p(0), to give 
j c = 2g 2 L exp(— \g\L) as required. At the critical current, we have that p(0) = 
2\g\Lexp(-\g\L) < 1, a = -2\g\Lexp(-\g\L) < 1, and \/3\ = 2g 2 > 1. 



2. The subcase \g\L < 1 

In the limit of \g\L -C 1, we have that p(x) is again close to 1 for all values of x. 
Consider the limiting expressions, eqn. (|T5|), for the parameters a and (3: we see that taking 
a — > j3 ~ p — p(0)/2 and hence j 2 — > (p — p/2) 2 p(0) leads to a local maximum for the 
supercurrent as the solutions loses its applicability above such values of j. In this limit, we 
have k ^> 1 and the form (|26|) follows for the density. The normalization condition, (|17|), 
gives 

(H + i + ^(o))§ = (H + |p(o))§. (C2) 

where we have used |p|£ 2 1, which we will find to hold given \g\L 1. We see that the 
normalization condition is satisfied automatically as p(0) ~ 1. We also have that p(L/2) ~ 
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p(0) + j 2 L 2 /4, and it may be checked retrospectively that p(L/2) ~ 1 for \g\L <C 1. The 
jump condition, (fl9|) , then leads to 

Jc ~ jfZ 2 " ~ 2jc ' 
and hence j c ~ (2\g\/ L) 1 / 2 as required. 



APPENDIX D: ATTRACTIVE IMPURITIES AT L 



2\G\ 



. In this section, we consider the special value of the separation, L = 2\g\, at which 
the critical current vanishes entirely in the attractive case. More specifically, we take L = 
2|g|(l + e) and consider the limit of very small e (which may be positive or negative). Using 
the limiting form of the density, eqn. (^, in the main text, we see that the boson-number 
fixing condition, eqn. ([T7|), leads to 



p(L/2) = p(0) + (1 - 0) 



,L 2 



1 - 



7-/T 



4 V" P(0) J ' 

Inserting this value for p(L/2) into the jump condition, eqn. QT$D, we find 

1 2 



J' 



p(0)p(L/2f 



-6 + P 







r p(o). 


+ 7 



+ 



p(0) p(L/2) 



At this point we are justified in substituting the values of p(0) and p(L/2) by their values at 
7 = 0: that is, we take p(0) as the solution of eqn. (p9|), while p(L/2) = p(0) + L 2 /4. This 
leads to the following expression for j 2 : 



J 



P(0) 



(2p(0) 



ep(0) 



1 



2p(0) \ 
p{L/2)) 



7 



ep(0)- 4p(0)-l 



2^7 
p{L/2))\ 



It may be verified that for e = 0, ie. L = 2\g\, the solution j = exists for 7 = 0, while 
as soon as 7 becomes non-zero, the solution for j becomes imaginary: the flow equations 
are unable to admit a well-behaved solution for any non-zero supercurrent. In contrast, for 
non-zero e, we may maximise j straightforwardly with respect to 7 to arrive at 



e 2 p(0) ; 



(2p(0)/p(L/2) + l)(4p(0) + 2p(0)/p(L/2) - 1) ' 
In particular, we see that j c vanishes linearly as |e| as e — > 0. 



(Dl) 
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